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Abstract 

The  long-term  goals  of  this  research  are  to  better  understand  and  accurately  model  low- 
to  mid-frequency  reverberation  in  shallow  water  environment.  Specific  goals  are  to 
develop  a  model  of  reverberation  for  conditions  (1-10  kHz,  ~20  m  water  depth,  ~  10  km 
range)  corresponding  to  the  ONR  Target  and  Reverberation  Experiment  performed  on  the 
spring  2013  (TREX2013),  develop  a  code  and  conduct  computer  simulations  with 
environmental  inputs  typical  for  the  chosen  location,  and  apply  this  model  to  analysis  of 
available  TREX2013  data.  This  report  presents  a  Green’s  function  modeling  approach 
that  allows  fast  estimations  of  volume  reverberation  in  complex  shallow  water 
environments.  A  simplified  first-order  version  of  the  approach  is  considered  to  show  how 
far-field  scattering  solutions  obtained  for  free  space  can  be  incorporated  into 
reverberation  in  complicated  bounded,  range-dependent,  and  stratified  environments.  A 
higher-order  modification  of  this  approach  is  considered  as  well,  using  a  MESB  (Multiple 
Eorward  Single  Backscatter)  approximation.  Application  to  analysis  of  TREX2013 
reverberation  and  tentative  model-data  comparisons  are  presented. 


Keywords 


Underwater  acoustic  propagation  and  reverberation,  Green’s  function  in  a  shallow 
water  waveguide,  volume  and  roughness  scattering,  data  inversion,  TREX13  data 
analysis  and  model-data  comparisons. 
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Objectives 


The  long-term  goals  of  this  research  are  to  better  understand  and  accurately  model 
low-  to  mid-frequency  propagation  and  reverberation  in  shallow  water  environment.  This 
would  result  in  development  of  a  robust  model  (potentially  a  tool)  for  prediction  and 
interpretation  of  reverberation  in  complex  (range-dependent,  potentially  3D-dependent) 
waveguides.  This  report  presents  a  Green’s  function  modeling  approach  that  allows  fast 
estimations  of  propagation  and  reverberation  in  shallow  water.  The  approach  also 
addresses  a  basic  science  question:  how  scattering  solutions  obtained  for  free  space  (or 
measured  in  direct-path,  short-range  conditions)  can  be  incorporated  into  long-range 
reverberation  in  bounded  and  complicated  environments,  such  as  shallow  water  or  deep¬ 
water  waveguides,  where  effects  of  multiple  interactions  with  boundaries  and  volume 
heterogeneities  must  be  taken  into  account,  providing  both  sufficient  accuracy  and  speed 
of  calculations. 

Specific  goals  of  this  research  are  to  develop  a  model  of  reverberation  for  conditions 
(1  -  10  kHz,  ~  20  m  water  depth,  ~  10  km  range)  corresponding  to  the  ONR  Target  and 
Reverberation  Experiment  performed  on  the  spring  2013  (TREX2013),  develop  a  code 
and  conduct  computer  simulations  with  environmental  inputs  typical  for  the  chosen 
location,  and  apply  this  model  to  analysis  of  available  TREX2013  data,  see  [TREX13 
web-site,  TREX13  Workshop  (2012,  2013,  2014)].  This  report  demonstrates  capabilities 
of  the  developed  codes  based  on  Green’s  function  modelling  approach  to  provide 
extremely  fast  predictions  of  reverberation  in  a  very  shallow  water  environment,  gives 
numerical  examples  for  several  relevant  scenarios,  discusses  possibilities  for  a 
quantitative  interpretation  of  reverberation  data,  presents  tentative  TREX2013  data- 
model  comparisons,  and  suggests  implications  for  future  research. 
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Project  Background 

A  proposal  for  this  project  was  written  in  response  to  the  ONR  call  for  suggestions  and 
comments  on  White  Paper  (2011)  “Mid- frequency  reverberation  measurements  with  full 
companion  environmental  support”,  which  defined  conditions  for  approaching 
reverberation  experiment  in  a  very  shallow  water  (~20  m  depth).  Scattering  mechanisms 
to  be  addressed  were  specified  as  well:  sea  surface  and  bottom  roughness,  heterogeneity 
in  water  column  (near  surface  bubbles,  fish  schools)  and  sediments  (e.g.  large  shells  and 
mud  inclusions).  White  Paper  emphasized,  however,  importance  to  address  applicability 
of  this  study  to  more  general  environments.  In  addition  to  the  general  problem  of 
reverberation,  specific  basic  science  issues  to  be  addressed  were  defined:  (1)  Mid¬ 
frequency  sound  scintillation  index  and  coherent  field,  (2)  Evolution  of  mid-frequency 
spatial  coherence  as  a  function  of  range  from  source,  (3)  Sediment  sound  speed  and 
attenuation  as  a  function  of  frequency,  (4)  Effects  of  spatial  variability  of  bottom 
scattering  strength,  (5)  Impact  of  sediment  ripple  fields  on  waveguide  propagation  and 
sediment  volume  scattering. 

The  proposal  for  this  project,  “Modeling  of  mid-frequency  reverberation  in  very  shallow 
water”,  submitted  in  2012,  can  be  briefly  described  as  follows.  The  main  goal  of  the  project 
was  to  enhance  the  modeling  component  of  the  approaching  experiment,  and  to 
contribute  to  development  of  non-traditional  approaches  to  shallow  water  reverberation, 
i.e.  to  certainly  support  the  idea  expressed  in  White  Paper  (2011),  that  a  new  area  of 
investigation  for  WPRM  (Wave  Propagation  in  Random  Media)  in  shallow  water  should 
be  suggested.  The  proposal  mentioned,  however,  that  for  such  very  shallow  water  and  to 
support  this  specific  experiment,  the  focus  of  such  investigation  should  be  broader  than  or 
even  different  from  the  one  in  the  mentioned  White  Paper,  where  effects  of  only  internal 
waves  were  emphasized  in  this  connection.  This  project  suggested  that  the  investigation 
most  relevant  to  very  shallow  water  environment,  where  sediment  is  a  critical  part  of  the 
propagation  channel,  should  be  focused  on  the  effect  of  3D  variability  of  the  sediment 
properties  along  with  corresponding  effects  of  propagation  and  3D  refraction  within  the 
seafloor. 

It  was  assumed  that  in  such  complicated  conditions  some  multiple  scattering  effects 
may  be  important,  and  suggested  to  address  these  effects  exploiting  an  approach 
developed  by  De  Wolf  (1971)  and  then  successfully  used  by  Ishimaru,  Tatarskii,  and 
others  for  description  of  electromagnetic  propagation  and  multiple  scattering  in  a 
turbulent  atmosphere.  Particularly,  the  approach  is  applicable  to  description  of  the  so- 
called  backscattering  enhancement  effect,  known  also  as  “weak  localization”  and 
“coherent  backscatter”  effects.  This  approach  is  well  recognized  in  WPRM  community, 
and  called  the  MESB -approximation  (Multiple -Eorward-Single-Backward-scatter).  While 
similar  terminology  appears  frequently  in  the  underwater  acoustics  literature,  such 
approach  for  conditions  of  shallow  water  reverberation  has  not  yet  been  adequately 
developed,  although  experimental  observations  exist,  Sabra  (2010). 

Eor  description  of  multiple-scattering  effects  in  acoustics  of  marine  sediments,  the 
MESB-approach  was  used  by  Ivakin  (1999),  and  then  it  was  reformulated  for  the  case  of 
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reverberation  in  a  shallow-water  waveguide  in  Ivakin  (2008).  The  main  complication  in 
this  case  (in  comparison  with  previous  work)  appears  because  of  the  multi-path 
conditions  for  propagation  in  a  waveguide.  Recent  developments  suggest  treating  this 
complication  by  direct  considering  both  the  heterogeneous  sediment  and  water  column  as 
parts  of  ocean  waveguide  Ivakin  (2011,  2012).  In  this  research,  the  approach  was 
proposed  to  be  developed  further  for  a  very  shallow  water  environment.  In  particular,  it 
was  suggested  to  perform  a  numerical  modeling  for  several  typical  scenarios  including 
both  water  column  and  seabed  as  possible  mechanisms  of  reverberation  to  better 
understand  their  relative  effects. 

One  of  practically  important  directions  of  this  research  is  development,  using  both 
analytical  and  numerical  modeling,  of  practical  relationships  between  reverberation  and 
environmental  parameters,  with  a  specific  goal  to  evaluate  the  sensitivity  of  measured 
reverberation  intensity  to  these  parameters,  to  determine  which  environmental 
parameters  or  mechanisms  of  scattering  are  critical  for  reverberation  in  chosen  acoustic 
and  environmental  very  shallow- water  scenarios.  This  may  help  to  reduce  the  number  of 
parameters  needed  to  be  measured  or  taken  into  account  to  provide  “full  environmental 
support”  of  the  field  experiment  proposed  in  White  Paper  (2011).  Another  practically 
important  application  of  this  modeling  would  be  evaluation  of  the  possibility  for  inferring 
environmental  parameters  from  reverberation  data.  New  techniques  and  algorithms  for 
such  inferring  were  expected  to  be  developed. 

This  project  was  proposed  for  a  two-year  period,  FY2013  and  FY2014.  During  the 
first-  year,  a  simplified  version  of  the  reverberation  model  was  planned  to  be  developed, 
with  a  focus  on  bottom  scattering  mechanisms,  to  facilitate  developing  codes  and 
conducting  pre-test  computer  simulations,  helping  with  planning  the  acoustic  experiment 
and  environmental  ground  truth  measurements.  Results  of  this  work  were  presented  in 
Ivakin  (2013),  and  Hefner  et  al  (2013a,  2013b).  During  the  second  year  of  work,  in 
FY2014,  a  more  general  model  of  reverberation  was  planned  to  be  developed  and  applied 
to  analysis  of  TREX2013  data  based  on  available  by  then  acoustic  and  environmental 
inputs.  Also,  it  was  suggested  to  consider  analytical  expressions  for  reverberation, 
methods  of  its  inversion  for  environmental  parameters,  and  applicability  of  results 
obtained  in  specific  TREX2013  conditions  to  more  general  shallow  water  scenarios. 
Some  results  of  this  work  are  presented  in  Ivakin  (2015a),  and  Hefner  et  al  (2015a, 
2015b). 
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Modeling  Approach 


This  report  presents  a  Green’s  function  modelling  approach  that  allows  fast 
estimations  of  reverberation  in  shallow  water.  The  approach  also  addresses  a  basic 
science  question:  how  far-field  scattering  solutions  obtained  for  free  space  can  be 
incorporated  into  reverberation  in  bounded  and  complicated  environments,  such  as 
shallow  water  or  deep-water  waveguides,  where  effects  of  multiple  scattering  and 
interaction  with  boundaries  and  volume  heterogeneities  must  be  taken  into  account. 

A  basic  approach,  its  simplified  first-order  version  and  initial  results,  are  based  on 
models  discussed  in  Ivakin  (2010,  2011,  2012,  2013).  A  higher-order  modification  of  this 
approach  exploits  ideas  of  a  MFSB  (Multiple  Forward  Single  Backscatter)  approximation 
developed  by  De  Wolf  (1971)  for  studying  electromagnetic  propagation  and  multiple 
scattering  effects  in  a  turbulent  atmosphere.  For  description  of  such  effects  in  acoustics  of 
marine  sediments,  the  MFSB-approach  was  applied  in  Ivakin  (1999),  and  reformulated 
for  the  case  of  reverberation  in  a  shallow-water  waveguide  in  Ivakin  (2008).  Its  further 
development  is  presented  in  Ivakin  (2015a). 


A  first-order  Green’s  function  approach 

Consider  the  total  field  (acoustic  pressure)  as  a  sum  of  the  unperturbed  field  and  a 
first-order  (single-scattered)  field  expressed  through  correspondent  components  of  the 
full-field  Green’s  function, 

(1) 

with  a  coefficient,  A,  introduced  so  that  Green’s  function  near  the  source  represents  a 
spherical  wave  of  the  unite  magnitude,  i.e.  |F-f"|G(f,f") — >1 .  Analytical  results 

and  numerical  solutions  for  zeroth-order  Green’s  function  are  currently  available  for 
rather  complicated  background  (or  unperturbed)  media,  such  as  range-dependent  shallow 
water  or  deep-water  waveguides,  and  are  based  on  various  propagation  models,  such  as 
PE,  normal  modes,  wave-number  integration,  or  various  ray-based  approximations. 

A  solution  for  the  single-scattered  field  is  known,  see  e.g.  Ivakin  (2008,  2010,  2011, 
2012,  2013),  which  then  results  in  an  integral  expression,  Ivakin  (2011,  2012) 

(/i(F,F'))  =  |A|"f(po(F)//7o(7))"M,(f)|Go(7,F)Go(f,F')fjV  = 

(2) 

=  |A|  ^jMy^^(r)/g(r,r')/o(r,r")d^r 

I  |2 

where  /oi  =  f’oi  the  zeroth-  (unperturbed)  and  first-order  (scattered)  intensities, 
Pair)  is  the  density  in  unperturbed  medium,  =  {p^fr' ) !  p^ff-ff  M  y ,  and  My  is  a 
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volume  scattering  coefficient,  or  the  scattering  cross-section  per  unit  volume.  Limitations 
of  Eq.(2)  result  from  assumption  that  the  spectrum  of  heterogeneity  is  a  smooth  enough 
function,  allowing  neglecting  bistatic  scattering  effects  by  ignoring  difference  in  wave 
vectors  appearing  in  the  stratified  environment  with  both  down-going  and  up-going 
waves,  Ivakin  (2011,  2012).  Despite  of  its  simplicity,  Eq.(2)  is  rather  general  and  can  be 
applied  to  scattering  near  marine  boundaries,  seabed  or  sea  surface,  or  near  an  arbitrary 
reference  interface  within  marine  environment,  using  a  general  procedure  described  in 
Ivakin  (2012). 

A  specific  expression  for  My  can  be  obtained  from  a  far-field  solution  for  scattering 

in  free  space  with  a  specified  type  of  heterogeneity,  continuum  or  discrete.  Eor 
backscatter  from  continuum  heterogeneity  caused  by  local  spatial  fluctuations  of  the 
density  and  sound  speed,  we  have,  see  Ivakin  (2011,  2012), 

My=27ik^^^-{s^)  (3) 

where  -{pc)I{PqCq)-\  are  relative  fluctuations  of  the  impedance,  and  O  is  their 

power  spectrum  at  the  backscattering  Bragg’s  wave-vector.  In  the  case  of  discrete 
scatterers,  such  as  particles  and  objects  of  various  kinds,  a  general  expression  for  the 
incoherent  volume  scattering  strength  of  randomly  distributed  objects  is  of  the  form,  see 
Ivakm(2012), 


My=Cy{alv)  (4) 

where  cr  and  v  are  the  individual  scattering  cross  section  of  the  object  and  its  volume, 
their  ratio  is  averaged  (given  a  fixed  volume)  over  other  parameters,  such  as  size,  shape 
and  orientation,  and  Cy  «  1  is  the  total  volume  concentration  of  the  objects. 

Therefore,  the  approach  provides  a  simple  way  for  calculations  of  reverberation  as 
follows.  An  integral  expression  for  the  backscatter  intensity,  given  by  Eq.(2),  has  a 
factorized  integrand  comprised  of  two  kernels,  the  two-way  propagator  and  the  scattering 
kernel.  The  propagator  is  a  product  of  two  local  intensities,  each  being  defined  along  one 
of  the  two  ways  of  propagation.  Eor  rather  complicated  environments,  generally  depth- 
range-dependent,  the  local  intensity  can  be  calculated  using  available  models  and  codes, 
e.g.  PE,  normal  modes,  ray-based,  or  other  approximations. 

The  scattering  kernel  is  defined  as  a  local  volume  scattering  coefficient  and  exploits 
simple  first-order  solutions  for  far-field  scattering  from  a  heterogeneous  volume  in  free 
(unbounded)  space.  It  can  be  specified  for  sea-water  column  and  seabed  with  continuum 
heterogeneity,  such  as  spatial  fluctuations  of  density  and  sound  speed,  and/or  discrete 
randomly  and  sparsely  distributed  targets,  such  as  gas  bubbles,  fish,  shells,  lenses-like 
inclusions,  oil  droplets,  solid  hydrate  particles,  and  others,  Ivakin  (2011,  2012). 

Moreover,  the  scattering  kernel  can  include  a  component  due  to  contribution  of 
roughness  at  an  arbirary  number  of  interfaces  at  depths  z  =  Zj  as  follows 

My^  (f )  =  m^^My  +  -  Zj)  (5) 

j 
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(6) 


(V) 


where  m  =  p !  and  n  =  c^l  c  are  relative  density  and  refraction  index  of  the 

background  medium,  M and  Dj  are  the  roughness  scattering  coefficient,  the 

power  spectrum,  and  a  contrast  factor  at  j-th  interface,  respectively.  In  more  detail,  this 
case  is  described  in  Ivakin  (1998,  2015c). 

Note,  that  this  approach  provides  estimations  of  bottom  reverberation  without 
calculations  of  the  equivalent  surface  scattering  strength  (although  may  include  it  as  a 
particular  case,  Ivakin  (2012,  2013,  2015a).  The  approach  is  applicable,  generally,  for  an 
arbitrary  distribution  of  scatterers  defined  by  their  volume  scattering  coefficient  My ,  and 

roughness  scattering  coefficient  ,  see  Eqs.(3-5).  Therefore,  it  allows  a  fast  estimation 
of  potential  contributions  of  different  mechanisms  of  scattering  with  arbitrary  strengths 
and  locations. 

A  modified  higher-order  approach 

Integral  expression  (2)  for  the  reverberation  intensity  can  be  generalized  by  treating 
the  propagating  kernel  stochastically  and  accounting  for  multiple  forward-scatter  effects, 
using  MFSB  approximation  developed  by  De  Wolf  (1971).  According  to  this  approach, 
the  total  field  is  presented  as  a  sum  of  a  multiple  forward-scattered  field  and  a  single- 
backscattered  field,  or  through  correspondent  components  of  the  full  Green’s  function. 


(8) 


Then  one  obtaines  the  integral  expression 


(9) 


1^1  ^  j  {^Vef  h  )h  ^2 


2 


where  I ^  ^  “  If*/  *1  multiple  forward-  and  single  back-scattered  intensities. 

Despite  of  apparent  similarity  to  the  first-order  expressions  (1,2),  Eqs.(8,9)  represent  a 
substantial  modification  and  improvement.  The  modified  approach  allows  taking  into 
account  such  multiple  scattering  effects,  as  the  scintillations  of  propagated  intensity,  the 
two-way  propagation  coherence,  and  the  related  “backscattering  enhancement”  known 
also  as  “coherent  backscatter”  and  “weak  localization”. 
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To  show  how  and  where  this  improvement  appears  and  ean  be  used,  reeall  an 
important  energy  eonserving  advantage  of  the  forward  seattering  propagator,  resulting  in 

{l  ~  Iq  .  Then  for  the  baekseatter  from  uniform  seatterer  distributions  one  obtains 


(4)  =  (4(4h))  ~  H  ^\My^(r)I^(r,r^)[l  +  K,(r,ri)]d^r 


(10) 


where  Kj  is  the  scintillation  index.  The  effect  of  backscattering  enhancement  appears 
here  automatically  because  Kj  >0. 

Consider  now  applicability  of  this  modified  approach  to  analysis  of  the  seabed 
reverberation  in  complicated  conditions,  such  as  in  the  TREX13  environment,  where 
marine  sediments  have  a  substantial  lateral  variability.  To  apply  Eq.(9),  a  general 
procedure  can  be  used  similar  to  described  in  Ivakin  (2012).  Consider  spatial  variables  of 
integration  in  (9)  of  the  form  r  -  (r^,^),  iy  -  >  and  introduce  normalizing  functions, 

for  example  as  incident  forward  propagation  wave  functions  on  the  bottom  surface, 
Go,  (4,1^2),  so  thatG^(r,rj  2)  ^  GQf(r^,r^  2)gi  2(^) .  Then  Eq.(9)  can  be  presented  in  the 
form 


(4  (4  C ))  =  l^r  J  (^5(4  ))(|Go  /  (4 ,  q  )Go^  (iy ,  4  )|"  y  V  (11) 

where  Mg  is  a  local  equivalent  interface  scattering  cross  section,  or  dimensionless 
bottom  scattering  strength, 


(12) 


Eqs.(9-12)  can  be  used  for  analysis  of  volume  reverberation  in  a  complex  TREX13 
environment  taking  into  account  heterogeneity  of  both  water  column  and  the  seabed.  Eor 
instance,  TREX13  sediments  were  very  heterogeneous,  and  had  a  substantial  lateral 
variability  of  the  sediment  type  and  composition,  mostly  due  to  variations  in  mud  and 
shell  content,  Hefner  and  Tang  (2014). 

The  observed  variability  of  sediment  composition  (e.g.  change  from  well-sorted 
medium  sand  to  muddy  or  shelly  sand)  may  result  in  noticeable  variations  in  sound  speed 
and  attenuation.  Particularly  important  in  this  case  are  changes  in  local  critical  angle, 
affecting  local  reflection,  penetration,  and  scattering  properties.  Eor  instance,  alternations 
of  the  sediment  type  along  the  propagation  paths  may  have  a  strong  accumulation  effect 
on  the  propagator  kernel  in  (9-1 1),  resulting  in  the  propagation  intensity  scintillations  and 
corresponding  backscattering  enhancement,  see  Eq.(lO).  The  alternations  across  the  path 
result  in  an  effective  averaging  of  the  scattering  kernel  in  (11)  given  by  the  local  bottom 
scattering  strength  Mg  ,  which  is  very  sensitive  to  changes  of  the  sediment  type, 
especially  at  shallow  grazing  angles  (below  the  critical  angle).  Therefore,  such  lateral 
variability  may  significantly  affect  the  reverberation  at  long  ranges.  A  quantitative 
analysis  of  the  effect  based  on  this  approach  is  being  planned  as  a  part  of  future  work. 
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Numerical  Examples 


Examples  here  are  to  demonstrate  capabilities  of  the  first-order  Green’s  function 
modeling  approach  to  provide  fast  estimations  of  reverberation,  to  show  three  critical 
steps  in  the  algorithm  of  these  estimations,  and  to  help  analyze  potential  contributions  of 
different  mechanisms  of  scattering  in  shallow  water  environments  (heterogeneous  water 
column,  rough  and  heterogeneous  bottom).  Several  scenarios  with  different  types  of  a 
layered  mud/sand  bottom  are  considered,  which  are  applicable  to  very  shallow  water,  as 
well  as  to  more  general  shallow  water  environments,  see  e.g.  ONR  Workshop  (2015). 
Input  parameters  for  bottom  are  chosen  to  be  representative  and  typical  for  mud  and  sand 
sediments. 


Environmental  scenarios  and  input  parameters 

In  the  following  numerical  examples,  Green’s  function  magnitude,  as  a  function  of 
range  and  depth  (specified  for  given  very  shallow  water  range-independent  environment), 
was  calculated  using  a  PE  code,  Tang  (2014),  with  input  parameters  as  follows: 

Erequency  -  3.5  kHz; 

Water  depth  -  20  m; 

Source/recever  depth  -19  m; 

Three  bottom  scenarios  were  chosen  to  represent  sand/mud  sediments  with  different 
thickness  of  the  mud  layer  (h)  placed  above  sand  halfspace: 

(1)  sand  halfspace  (h=0), 

(2)  0.5  m  of  mud  above  sand  halfspace, 

(3)  3.5  m  of  mud  above  sand  half  space. 

Calculations  for  the  penetration  field  within  the  bottom  were  made  for  depths  down  to 
5  m  below  bottom  surface. 

Acoustic  parameters  for  water,  mud,  and  sand  are  taken  as  follows: 


water 

mud 

sand 

Sound  speed  [m/s] 

1530 

1500 

1670 

Density  [g/cm^] 

1.0 

1.6 

2.0 

Attenuation  (loss  tangent),  S 

0 

0.002 

0.01 
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Numerical  simulations  in  steps 

The  algorithm  of  fast  estimations  of  shallow  water  reverberation  ineludes  three  critical 
steps.  The  first  step,  according  to  Eq.(2),  is  to  evaluate  the  propagation  kernel  as  a 
function  of  locations  of  the  source/receiver  and  scattering  point.  For  this,  the  Green’s 
function  magnitude  should  be  pre-calculated  in  the  waveguide  of  interest.  Importantly, 
the  algorithm  requires  calculations  only  at  one  (central)  frequency.  These  calculations  are 
illustrated  in  Figures  1-3. 


Figure  1.'  Intensity  of  3.5  kHz  PE-propagation  field  in  shallow  water  (of  20m  depths)  with 
sand  bottom.  The  source  is  1  m  above  bottom. 


PE  3D-corrected  primary  field  intensity  at  3.5  kHz,  0.5m  mtxl  over  sand 


Range  (m) 


Figure  2:  Intensity  of  3.5  kHz  PE-propagation  field  in  shallow  water  (20m  depth)  with  a 
two-layer,  0.5m  of  mud  over  sand,  bottom.  The  source  is  1  m  above  bottom. 
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PE  3D-corrected  primary  field  intensity  at  3.5  kHz,  3.5m  mud  over  sand 


Range  (m) 


Figure  3.'  Intensity  of  3.5  kHz  PE-propagation  field  in  shallow  water  (20m  depth)  with  a 
two-layer,  3.5m  of  mud  over  sand,  bottom.  The  source  is  1  m  above  bottom. 


The  second  step  of  the  algorithm  is  illustrated  in  Figures  4-6.  It  is  seen  from  Figure  4, 
that  the  range-dependence  of  the  propagation  intensity  on  the  bottom  surface,  which  is 
important  for  estimation  of  bottom  reverberation,  is  a  strongly  oscillating  function.  For 
further  analysis,  such  calculations  require  a  kind  of  smoothing.  For  a  particular  kind, 
results  are  shown  in  Figure  5,  where  the  smoothing  was  performed  for  a  fixed  range  span 
L  -  100  m.  A  simple  rationale  for  such  smoothing  may  result  from  possible  uncertainty 
in  ranges/depths,  which  can  be  substantial  in  some  measurement  techniques,  e.g.  in  the 
case  where  source  or/and  receiver  locations  are  not  fixed.  For  instance,  such  technique 
was  used  in  the  TREX  propagation  measurements  with  a  fixed  vertical  array  in  mid-water 
(at  9-15  m  depths)  and  a  CW-source  moving  (i.e.  not  fixed)  at  a  10  m  depth,  so  that  a 
range  dependence  of  the  Transmission  Loss  (TL)  was  measured.  It  was  noticed  that,  for 
regions  with  sand  bottom  covered  by  a  mud  layer  (about  20cm  thick)  and  regions  with 
sand  bottom  (free  of  mud),  TL  is  essentially  the  same.  This  effect  is  easy  to  see  in  Figure 
6,  where  even  a  50  cm  mud  layer  is  shown  to  be  not  sufficient  for  a  noticable  effect  on 
propagation  intensity  averaged  over  100  m  ranges  and  9-15  m  depths  spans.  Another 
reason  for  smoothing  over  a  range  span  will  be  discussed  later  in  connection  with  TREX 
reverberation  data  analysis. 
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Figure  4;  Range-dependence  of  Transmission  Loss  at  the  bottom  surface  (at  20m  depth), 
3.5  kHz  PE-propagation  field  for  three  different  types  of  bottom:  sand  half- space 
(red),  0.5m  of  mud  over  sand  (blue),  and  3.5m  of  mud  over  sand  half- space 
(green).  The  source  is  1  m  above  bottom.  Also  shown  are  and  R^  range- 
dependences. 


TL  at  depth  20m,  averaged  over  a  range  span  100  m,  scenarios  1-3 


-  --  1 

f 

-  --  - 

/ 

—  - 

-  — 

■ 

■ 

- 0.5m  mud  over  sand 

- 3.5m  mud  over  sand 

- R'^ 

1 000  2000  3000  4000  5000  6000  7000 

Range  (m) 


Figure  5.‘  Smoothed  over  100m  range  span  range-dependence  of  Transmission  Loss  at 
the  bottom  surface  (at  20m  depth),  3.5  kHz  PE-propagation  field  for  three 
different  types  of  bottom:  sand  half- space  (red),  0.5m  of  mud  over  sand  (blue), 
and  3.5m  of  mud  over  sand  half-space  (green).  The  source  is  1  m  above  bottom. 
Also  shown  are  R^'^  and  R^  range-dependences. 
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Range  (m) 


Figure  6:  Range-dependence  of  Transmission  Loss  in  mid-water,  averaged  (smoothed) 
over  100m  range  span  and  over  9- 15m  depths,  3.5  kHz  PE-propagation  field  in 
shallow  water  (20m  depth)  for  three  bottom  types:  sand  half-space  (red),  0.5m  of 
mud  over  sand  (blue),  and  3.5m  of  mud  over  sand  half- space  (green).  The  source 
is  1  m  above  bottom.  Also  shown  are  and  R^  range-dependences. 


The  third  step  of  the  algorithm  provides  extremely  fast  evaluation  of  reverberation  in  a 
shallow  water  waveguide  using  pre-calculated  for  this  waveguide  range-depth 
dependence  of  single-frequency  Green’s  function  intensity.  The  effect  of  the  frequency 
bandwidth  is  considered  at  this  step  as  follows.  According  to  Eq.(2),  the  integration 
should  be  performed  over  a  horizontal  cross-section  of  the  scattering  (ensonified) 
volume,  which  assumes  a  span  of  azimuthal  angles,  Acp  ,  (to  account  for  possible 
azimuthal  directivity,  e.g.  for  horizontal  arrays,  such  as  FORA).  Also,  the  integration 
over  scattering  area  presumes  another  span  of  ranges,  Ar,  which  may  be  significantly 
different  from  above  mentioned  smoothing  scale  L.  The  span  Ar  is  related  to  a  span  of 
double  (two-way)  propagation  time,  and  can  defined  (although  only  approximately)  in  the 
same  way  as  in  a  direct-path  (short-range)  scenario,  through  the  frequency  bandwidth, 
A/  ,  or  through  the  pulse  duration,  r  ,  as  follows:  Ar  =  cr  !  2  =  c !  Af  12 .  In  the  following 
calculations,  input  parameters  needed  at  this  step  are  as  follows:  Af  =  100 Hz,  Acp  =  2.6° 

.  This  gives  an  estimate  Ar  =  8  m. 

At  this  step,  reverberation  caused  by  bottom  volume  heterogeneity  and  rough 
interfaces,  as  well  as  volume  reverberation  from  scatterers  in  water  column,  can  be 
calculated,  generally,  for  an  arbitrary  distribution  of  scatterers  defined  by  their  volume 
scattering  coefficient  My  ,  and  roughness  scattering  coefficient  M ^  ,  see  Eqs.(3-7). 

Therefore,  this  allows  fast  estimating  and  comparing  relative  contributions  of  different 
scattering  mechanisms,  with  different  nature,  strengths,  and  locations. 
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Results  for  chosen  bottom  scenarios  are  presented  in  Figures  7-11,  showing 
reverberation  caused  by  different  types  of  scatterers:  heterogeneities  of  water  column 
(Figure  7),  bottom  volume  (Figures  8,  9),  and  bottom  interface  roughness  (Figures  10, 
11),  using  representative  values  for  input  parametrs.  My  and  . 


Water  Columr  Volume  Reverb  at  3.5  kHz,  M  =  1e-06  m'f  scenarios  1-3 


- sand,  no  mud  j 

- 0.5m  mud  over  sand 

- 3.5m  mud  over  sand  1 

1 

1 

_ 

1 000  2000  3000  4000  5000  6000  7000 

Range  (m) 


Figure  1:  Normalized  intensity  of  reverberation,  RL-SL,  at  3.5  kHz,  caused  by  volume 
scattering  in  water  column  with  My  —  10”^  m~* .  Shallow  water  (20m)  with  three 

different  types  of  bottom:  sand  half-space  (red),  0.5m  of  mud  over  sand  (blue),  and 
3.5m  of  mud  over  sand  half-space  (green).  The  source  is  1  m  above  bottom. 


Bottom  Volume  Reverb  at  3.5  kHz,  M  =  0.001  m”',  scenarios  1-3 


- sand,  no  mud  I 

- 0.5m  mud  over  sand 

- 3.5m  mud  over  sand  I 

' 

I 

1 000  2000  3000  4000  5000  6000  7000 

Range  (m) 


Figure  8.'  Normalized  intensity  of  reverberation,  RL-SL,  at  3.5  kHz,  caused  by  sediment 
volume  scattering  with  My  - 10”^  m  .  Shallow  water  (20m)  with  three  different 

types  of  bottom:  sand  (red),  0.5m  of  mud  over  sand  (blue),  and  3.5m  of  mud  over 
sand  (green).  The  source  is  1  m  above  bottom. 
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M  =  0.001  nr',  P  =  2,  c  =  1.0915,  6  =  0.01,  f  =  3.5  kHz,  h  =  10  m 


1 - Theory 

1  Urrtoert,  -27  dB 

1 

10“  10'  10^ 

Grazing  Angle,  deg 


Figure  9:  Angular  dependence  of  Bottom  Scattering  Strength  (Normalized  Short-Range 
Reverberation  Intensity),  RL-SL,  for  3.5  kHz,  sand  bottom  with  the  same  (typical) 
scattering  strength  per  unit  sediment  volume  used  in  Figure  8,  My  =  10~^  m  . 
Also  shown  (in  green)  is  Lambert  Law  at  -27  dB. 


Rougness  Reverb  at  3.5  W-te,  kY  =  1.3106e-05,  scenarios  1-3 


1 

- sand,  no  mud 

- 0.5m  mud  over  sand 

- 3.5m  mud  over  sand 

. 

- • - 

1  , 

1 000  2000  3000  4000  5000  6000  7000 

Range  (m) 


Figure  10.’  Normalized  intensity  of  reverberation,  RL-SL,  at  3.5  kHz,  caused  by  bottom 
interface  roughness,  with  Mj^  =1.3x10  Shallow  water  (20m  depth)  with  three 
different  types  of  bottom:  sand  half-space  (red),  0.5m  of  mud  over  sand  (blue),  and 
3.5m  of  mud  over  sand  half-space  (green).  The  source  is  1  m  above  bottom. 
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■  =  4.13.  w  =  0.00168m‘‘''2,  L  =  6m,  p  =  2.  c  =  1.0915,  8  =  0.00915,  f  =  3.5  kHz 


-  Theory  [ 

Lambert,  -27dB 

Grazing  Angle,  deg 


Figure  11.'  Angular  dependence  of  Bottom  Scattering  Strength  (Normalized  Short-Range 
Reverberation  Intensity),  RL-SL,  at  3.5  kHz,  sand  bottom  with  a  roughness 

scattering  coefficient  used  in  Figure  10,  M^,  =  1.3  x  10  ^ .  Also  shown  (in  green)  is 
Lambert  Law  at  -27  dB. 
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Application  to  TREX  Data  Analysis 


The  described  in  previous  sections  approach  was  refined  and  specified  for  TREX13 
environment,  geometry  and  other  input  parameters. 


TREX  environmental  scenario  and  input  parameters 

As  in  previous  examples,  Green’s  function  magnitude  was  calculated  using  a  PE  code, 
Tang  (2014),  but  for  a  slightly  different  configuration  and  parameters  corresponing  to  a 
Pekeris  waveguide  with  sand  bottom.  Parameters  for  these  calculations  were  chosen  as 
most  relevant  to  TREX13  conditions  and  available  by  now  acoustic  data,  Tang  (2015): 

Central  frequency  -  3450  kHz; 

Frequency  bandwidth  -  100  Hz; 

Water  depth  -19  m; 

Source  depth  -  17.8  m; 

Receiver  (FORA)  depth  -  16.9  m; 

Azimuthal  beamwidth  -  2.6  deg; 

Source  Fevel  (SF)  -  200  dB. 

Calculations  for  the  penetration  field  within  the  sand  bottom  were  made  for  depths 
down  to  2  m  below  bottom  surface. 

Acoustic  parameters  for  water  and  sand  are  taken  as  follows: 


water 

sand 

Sound  speed,  (m/s) 

1525 

1630 

Density  (g/cm^3) 

1 

2.0 

Attenuation  (loss  tangent),  5 

0 

0.01 
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TREX-specified  algorithm  for  numerical  simulations 

The  algorithm  for  calculations  of  reverberation,  specified  for  the  TREX  scenario,  also 
includes  three  steps,  but  with  some  distinctions  from  those  described  earlier.  The  first 
step,  in  contrast  with  previous  examples,  considers  a  more  general  case  of  bistatic 
geometry,  i.e.,  according  to  Eq.(2),  now  the  propagator  is  to  be  calculated  considering 
different  locations  of  the  source  and  receiver.  For  this,  the  Green’s  function  magnitude 
was  pre-calculated  for  refined  input  parameters,  given  above.  Again,  it  is  important  that 
the  algorithm  requires  calculations  only  at  one  (central)  frequency. 

The  second  step  of  the  algorithm  was  also  refined  taking  into  account  the  fixed-fixed 
source-receiver  geometry  of  the  TREX  13  reverberation  measurements  and  a  broad-band 
frequency  filter  (100  Hz)  chosen  for  the  TREX  13  reverberation  data  processing,  Yang  et 
al  (2014).  This  case  requires  a  technique  for  smoothing  strong  oscillations  of  the 
propagation  intensity,  which  differs  from  that  used  in  previous  examples.  It  takes  into 
account  a  range-dependent  scales  in  interferencial  structure  of  the  broad-band  field  in 
shallow  water  waveguides.  In  this  case  the  range  span  for  the  smoothing  is  range 
dependent  as  well,  and,  in  contrast  with  previous  numerical  examples,  is  defined  as 
follows:  L  =  rAf  /  /  .  More  details  on  this  technique  are  described  in  Ivakin  (2015b). 

The  third  step  of  the  algorithm  provides  extremely  fast  evaluation  of  TREX  13 
reverberation  using  range-depth  dependence  of  single- frequency  Green’s  function 
intensity  pre-calculated  for  this  environment  and  smoothed  in  correspondence  with  the 
above-described  technique.  The  integration,  according  to  Eq.(2),  was  performed  over  a 
horizontal  cross-section  of  the  scattering  (ensonified)  volume,  with  a  span  of  ranges,  Ar , 
and  a  span  of  azimuthal  angles,  Acp ,  defined  in  the  same  way  as  in  previous  examples, 
i.e.  A(p  =  2.6°  and  Ar  =  8m. 


TREX  Reverberation:  tentative  model/data  comparison 

Eor  further  analysis,  reverberation  caused  by  bottom  volume  heterogeneity,  as  well  as 
volume  reverberation  from  scatterers  in  water  column,  was  calculated  and  a  tentative 
model/data  comparison  for  a  sub-set  of  TREX  reverberation  data  was  performed.  The 
data  sub-set  used  here  is  limited  to  a  published  example  of  TREX  reverberation  data, 
Hefner  and  Tang  (2014),  presented  in  some  more  detal  in  Yang  et  al  (2014),  see 
TREX13  Web-site.  This  example  includes  only  one  frequency  sub-band  (/  =3400^3500 
Hz),  one  azimutal  look  angle,  cp  =  129°  ,  and  only  two  runs.  Run  #17  (23  April  2013)  and 
Run  #79  (9  May  2013).  The  data  are  shown  in  Figure  12,  where  a  normalized 
reverberation  intensity,  Reverberation  Level  (RE)  corrected  for  Source  Level  (SL),  i.e. 
RE  -  SL,  is  presented. 

It  is  seen  that  the  data  shown  in  Figure  12  are  somewhat  different.  This  difference 
might  represent  an  effect  of  environmental  changes  resulted  from  a  storm  occurred  in  the 
area  around  1  May  2013.  The  following  calculations  and  results  of  model-data 
comparisons,  presented  in  Figures  13-19,  give  a  tentative  interpretation  of  this  effect. 
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TREX  Re\erberat!on  3.4-3. 5  kHz,  averaged,  before  storm  (red),  run  17,  and  after  (blue),  run  79 
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Figure  \2:  TREX13  normalized  reverberation  intensity,  RL-SL,  at  3.4-3. 5  kHz,  for  two 
runs,  before  (red)  and  after  (blue)  storm. 


The  model  suggests  that  reverberation  is  caused  by  volume  scattering  either  in  water 
column  or  in  the  sediment.  The  contribution  of  bottom  roughness  can  be  considered  in  a 
way  similar  to  shown  in  previous  numerical  examples  and  is  not  included  in  this  report.  A 
separate  paper,  Ivakin  (2015c),  which  discusses  this  case  in  detail,  is  in  preparation. 
Typical  values  for  the  only  free  parameter  of  the  model.  My ,  were  used  to  provide  a 

reasonably  good  fit  between  model  results  and  the  data. 

First,  reverberation  from  water  column  was  considered.  The  results  for  before  storm 
model-data  comparison  is  shown  in  Figure  13.  Generally,  it  shows  that  scattering  in  water 
column  can  provide  a  good  fit  to  data  at  long  ranges  (more  than  6  km).  However,  to  fit 
the  data  level,  volume  scattering  strength  in  water  column  should  be  somewhat  higher 
than  typical.  Unfortunately,  there  is  no  any  available  TREX  ground  truth  data  for 
scattering  in  water  column  at  these  ranges. 

Next  case  to  consider  is  volume  scattering  in  the  sediment.  As  an  initial  step  on  this 
way,  the  bottom  volume  scattering  coefficient  was  assumed  to  be  range-independent.  The 
results  are  shown  in  Figure  14.  Then  a  possibility  of  a  smooth  range-dependence  of 
My(r)  was  considered  to  reduce  the  model-data  difference.  This  resulted  in 

estimated/inferred  range-dependence  My(r)  shown  in  Figure  15.  Corresponding 
corrected  reberberation  model-data  comparison  is  shown  in  Figure  16. 
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Water  Column  Reverb  at  M  =  3e-07  m  ',  TREX  data  3.4-3.5  kHz,  a  =  0.033,  Run  17,  averaged  over  13  pings,  (|)  =  2.2° 
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Figure  13.'  TREX13  normalized  reverberation  intensity,  RL-SL,  at  3.4-3. 5  kHz,  measured 
before  storm  (red),  compared  to  model  result  calculated  for  scattering  in  water 
column  with  volume  scattering  coefficient  My  =  3  x  10”^  m  . 


TREX  data  3.4-3.5  kHz,  Run  17,  averaged  over  13  pings,  i|i  =  2.2°,  Initial  Model,  =  0.002  m  ',  a=  0.033 
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Figure  \A:  TREX13  normalized  reverberation  intensity,  RL-SL,  at  3.4-3. 5  kHz,  measured 
before  storm  (red),  compared  to  that  calculated  for  volume  scattering  in  sand 
bottom  ( range -independent)  with  My  =  0.002  m  . 
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Inferred  Range-Dependent  Volume  Scattering  Coefficient 


Figure  IS;  A  correction  made  for  range-dependence  of  local  volume  scattering  strength 
in  TREX13  sand  bottom  before  storm. 


Corrected  Model  accounting  for  range-dependence  of  M  vs  TREX13  Data 
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Figure  16.'  TREX13  normalized  reverberation  intensity,  RL-SL,  at  3.4-3. 5  kHz,  measured 
before  storm,  compared  to  that  calculated  for  volume  scattering  in  sand  bottom 
with  corrected  range-dependent  volume  scattering  strength  shown  in  Eigure  15. 
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A  similar  procedure  was  performed  for  analysis  of  after  storm  reverberation.  The 
results  for  initial  model  data  comparison  (for  range  independent  My )  are  shown  in 

Figure  17.  It  shows  only  a  slight  discrepancy  at  small  ranges,  which  was  easily 
compensated  by  a  correspondingly  slight  range-dependence  of  sediment  volume 
scattering  coefficient  shown  in  Figure  18.  Then  model-data  comparison,  shown  in  Figure 
19,  becomes  reasonably  good. 


TREX  data  3.4-3.5  kHz,  Run  79.  averaged  over  18  pings,  (|i  =  2.2°,  Initial  Model,  M  =  0.02  m'',  a  =  0.033 
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Figure  17.'  TREX13  normalized  reverberation  intensity,  RL-SL,  at  3.4-3. 5  kHz,  measured 
after  storm  (blue),  compared  to  that  calculated  for  volume  scattering  from  sand 
bottom  ( range-independent),  )  with  My  =  0.02  m  . 


Inferred  Range-Dependent  Volume  Scattering  Coefficient 


Figure  18.'  A  correction  made  for  range-dependence  of  local  volume  scattering  strength 
in  TREX  13  sand  bottom  after  storm. 
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Corrected  Model  accounting  lor  range-dependence  of  M  «  TREX13  Data 
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Figure  19:  TREX13  normalized  reverberation  intensity,  RL-SL,  at  3.4-3. 5  kHz,  measured 
after  storm  (blue),  compared  to  that  calculated  for  volume  scattering  in  sand 
bottom  with  corrected  range -dependent  local  volume  scattering  strength  shown  in 
Figure  18. 


Comparison  of  Figure  15  and  Figure  18  may  help  to  understand  nature  of  the 
environmental  impact  of  the  storm  in  the  area  that  significantly  affected  reverberation. 
Such  impact  could  result  from  dynamical  mixing  of  the  sediment  and  consequent 
smoothing  of  the  range  dependence  of  the  strength  of  surficial  sediment  heterogeneity 
along  the  observation  path.  Importantly,  same  conclusions  would  be  reached  if  the  range- 
dependence  of  local  roughness  scattering  is  considered  using  Eq.(5).  In  more  detail,  this 
will  be  described  in  a  separate  paper,  Ivakin  (2015c). 

In  addition,  to  ensure  that  the  two  runs,  #17  and  #79,  see  Figure  12,  used  here  for 
model-data  comparisons,  are  indeed  representative  of  a  more  comrehensive  data  set, 
reverberation  data  from  five  more  runs  which  recently  became  available,  Tang  (2015), 
were  considered.  They  are  added  to  previous  two  runs  and  shown  in  Figure  20  and 
confirm  the  general  effect  of  the  storm  seen  in  Figure  12.  This  suggests  that  above-shown 
model-data  comparison,  as  well  as  conclusions  based  on  this  comparison,  will  stay 
reasonable  for  a  more  comprehensive  TREX13  data  analysis. 


23 


Distribution  Statement  A:  Approved  for  public  release;  distribution  unlimited. 


Figure  20/  TREX13  normalized  reverberation  intensity,  RL-SL,  at  3.4-3. 5  kHz,  measured 
before  (red)  and  after  (blue)  storm.  Five  more  runs  added  to  Figure  12,  showing 
that  previously  used  two  runs  are  typical. 
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Concluding  Remarks  and  Implications  for  Future  Research 


This  report  has  presented  a  modeling  approach  that  allows  fast  estimations  of  volume 
reverberation  in  complex  shallow  water  environments.  Computer  simulations  with 
environmental  inputs  typical  for  shallow  water  have  been  conducted.  Capabilities  of  this 
modelling  approach  and  developed  codes  to  provide  a  reasonable  interpretation  of 
TREX2013  reverberation  data  are  demonstrated  and  discussed. 

Based  on  these  results,  future  research  may  address  several  science  issues,  which 
appeared  to  be  important  for  better  understanding  TREX2013  data: 

•  To  account  for  sediment  lateral  variability,  mud  patches/strips,  using  more 
advanced  PE  codes  for  Green’s  function  (here  only  the  case  of  Pekeris  waveguide 
was  considered); 

•  To  analyze  environmental  and  acoustic  data  to  get  inputs  describing  variability  of 
shell  content,  sub-bottom  interface  spectra,  heterogeneity  in  mud; 

•  To  extensively  support  TREX13  data  analysis,  including  data  obtained  for  whole 
mid- frequency  range,  1-10  kHz,  and  all  other  azimuthal  angles  (here  only  a  subset 
of  data  at  3. 4-3. 5  kHz  and  one  azimuthal  direction  was  considered). 

Several  journal  papers  are  planned  to  be  prepared  for  submission  in  near  future,  based 
on  results  and  ideas  of  this  research,  briefly  discussed  and  presented  in  this  report.  One  of 
them,  Ivakin  (2015b),  soon  to  be  submitted,  includes  results  shown  in  Eigures  12-19  of 
this  report  and  use  data  currently  available  on  TREX2013  web-site  (Eigure  12).  More 
papers  are  in  preparation,  tentatively  shown  as  Ivakin  (2015c,  2015d,  2015e). 
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